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Abstract: 

Neutrinos are one of the major puzzles in modern physics. Despite measurements of mass differences, the 
Standard Model of particle physics describes them as exactly massless. Additionally, recent measurements from 
both particle physics experiments and cosmology indicate the existence of more than the three Standard Model 
species. Here we review the cosmological evidence and its possible interpretations. 
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1 Introduction 

Neutrinos are the lightest of the known massive particles and one of the least understood. The cosmic neutrino 
background has a number density of 112cm -3 (per specie), which means that billions of neutrinos pass through 
our bodies every second, so they are by no means rare. Yet we do not know their basic properties, such as their 
mass and how many species there are. 

1.1 Particle physics neutrinos 

The discovery of th e Higgs boson at the Large Hadron Collider (LHC) consolidates the power of the Stan dard Model 
of particle physics (| ATLAS Collaboration Aad et al]|2012l : [c"MS Collaboration Chatrchvan et al.ll2012l '). However, 
there are still unsolved issues in the neutrino sector, with their masses remaining unexplained and uncertainty 
about the actual number of species. 

The Standard Model contains three species of mass-less neutrinos that interact only through the weak inter- 
action. By measuring the decay width of the Z° boson the Larg e Electron Positron C ollider (LEP) measured the 
number of weakly interacting neutrinos to be N v = 2.92 ± 0.05 (Bcringcr et al. 20l2), corresponding to the three 

generations of leptons namely the electron, the muon and the tau pa rticles. 

Observations of neutrino disappearance, e. g. from v e sou rces ((Cleveland et al ] 119981 : iFukuda et aD Il998l ). 



and neutrino appearances e.g. in beams ( Abe et all 12011). have lead to the framework of neutrino mix- 



ing/oscillations, where the interaction eigenstates {v e .,v li ,v T ) are described as superpositions of mass eigen- 
states (v\, v^, vz)- Global fits to data from experiments using solar, atmospheric, and reactor neutrinos have 
measured mass differences between the three species to be Am| 2 = |(2.43jio'jjj| ) x 10 _3 | eV 2 and Am| a = 
(7.50 ± 0.20) x 10" 5 eV 2 l|Beringer et all I2012T I thus requiring at least two mass eigenstates to have non-zero 
masses. 

From studying the end-point of the electron energy distribution in /3-decays, the Mainz and Troitsk experiments 
have limited the electron neutrino to hav e a mass smaller than 2.3 eV and 2.05 eV respectively (95 % confidence 
level E et al.1 120051 : II et al.ll201 If ). For the specific case of neutrinos being Majorana particles (i.e. their 
own anti-particles), the Heidelberg- Moscow experiment has limited the effective neutrino mass to be less than 
0.35 eV (90% confidence level) using neutrino-less double /3-decay (Klapdor-Klei ngrothaus fc Krivosheinal [20061 ') 
with similar results from EXO (lAuger et al.ll2012T ). However, no present-day experiment has the sensitivity to 
measure the absolute neutrino mass. The upcoming KArlsruhe TRItium Neutrino (KATRIN) experiment is 
expected to improve the sensitivity to 0.2 eV (90% confidence level) from /3-spectroscopy, but this upper limit 
remains far above the observed mass differences. 

The sign of Am32 is unknown, which means that we cannot know whether the ordering of the mass states 
is mi < 77i2 << 7713, similar to the leptons, where m e < m M << m T . This is called normal hierarchy. If 
7773 << 777i < m2, the scenario is called inverted hierarchy. 

For the short baseline (SBL) oscillations there are some tensions between the various experiments allowing 
for, or even preferring, the existence of additional neutrino species that do not interact through the weak force, 
and are consequently denoted sterile. Depending on the exact analysis, the preferred scenario includes either one 
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or two sterile neutrinos in addit i on to the three normal ones (3+ 1 and 3+2), with sterile neutrino masses on the 
order of 1 eV (|Kopp et alj|201ll ; iMention et aL|[201ll ; lHuber|[201ll ; iGiunti fc Lavededl20Tlh . 



1.2 Cosmological neutrinos 

Observations of the Cosmic Microwave Background (CMB) have, since the 5-year data release from the Wilkinson 
Microwave Anisotropy Probe (WMAP), displayed a mild preference for an excess of relativistic energy density 
at the time of decoup ling on top of what is provided by the photon density and three Standard Model neutrinos 
(|Komatsu et al.ll2009l ). This extra radiation has been dubbed 'dark radiation' and is often parametrised in terms 
of extra neutrino species such that N e g = N v + AiV. The most striking result is from the Atacama Cosmology 
Telescope (ACT) of N cH = 5.3 ± 1.3 (|Dunklev et al.ll201ll ) which is shown in Fig. [T] together with a selection of 
cosmological analyses from the literature. There is a clear trend favouring N a ff> 3 across various models and data 
sets. 

./Veff is often referred to as the effective number of neutrino species or the number of relativistic species, 
but basically it is a parametrisation of anything that would change the expansion rate at the time of decoupling. 
Consequently fractional values of N B g do not actu ally require a non-integer number of neutrino species. Possibilitie s 
include other exotic particles such as axions (e.g. [Melchiorri e"t~a l. 2007; Nakavama et al.ll201ll:]Erken^t_alJl2Qld), 
decay ing particles or fields (jlchikawa et alj 120071 ; iNakavMn^etalT 20111; lFischler*fc*^eversl |2011| ; iBoefrrnet al.l 
120121) . gravity waves dSmith et al.ll2006l). extra dim ensions l|Flambaum fc Shurvakll2006l ) "and early dark energy 



( Calabrese et al.ll201ll ; lGagnon &i Leseourgues!l201ll . and references therein). 



It is tempting to connect the cosmological preference for N c g> 3 with the possible additional species from 
oscillation experiments. In this review we focus on the physical effects (Sec. [2]), the cosmological observations 
(Sec. the possible origin of these effects in terms of particle physics and other explanations (Sec. 2J), and what 
to expect in the future (Sec.JS]). 

Th roughout the paper we use a fiducial cosmology with the best fit WMAP7 parameters from lKomatsu et al.1 
(120111 ) for illustrations: Q,b,o = 0.0451 (baryon density today), fi Ci o = 0.226 (matter density today), £Ia,o = 0.729 
(dark energy density), Ho = 70.3 km s _1 Mpc -1 (Hubble parameter today), n B — 0.966 (index describing the tilt 
of power spectrum of primordial density fluctuations), A^(fc = 0.002 Mpc -1 ) = 2.42 x 10 -9 (amplitude of the 
primordial power spectrum at k — 0.002 Mpc" 1 ), and unless otherwise stated w = —1 (equation of state of the 
dark energy). All quoted uncertainties on N e s are 68% confidence levels, and all upper limits on ^tn„ are 95%. 



2 Early Universe effects 

In order to understand the observable effects of dark radiation parametrised as N e ff, it is illustrative to follow the 
normal neutrinos during the evolution of the early Universe depicted for the fiducial cosmology in the left column 
of the timeline in Fig. [2] 



Neutrino decoupling 

(T„ dcc » 1.1 MeV, z udcc « 4.5 x 10 9 ) 

Neutrinos only interact through the weak interaction and consequently when the early Universe temperature drops 
below 1 MeV, the interaction rate per neutrino drops below the expansion rate, and the neutrinos decouple from 
the thermal plasma. For neutrino masses below 1 eV the neutrinos are still ultra-relativistic when they decouple 
and they retain a relativistic Fermi-Dirac phase-space distribution. Their energy distribution is described by a 
blackbody spectrum with a single temperature, and they may have a non-zero chemical potential (related to lepton 
asymmetry) . 



Big Bang Nucleosynthesis 
(T nuc » 0.8 MeV, z nuc » 2.6 x 10 9 ) 

The formation of deuterium (D) and helium (He) nuclei is affected by the neutrinos in two ways. Firstly, electron 
neutrinos participate in the charged current weak interactions that determine the neutron to proton ratio (e.g. 
!/ e +n-H>p + e~), and secondly, they increase the expansion rate as 

H 2 (z)~^( Pj + p„), (1) 

where H(z) is the Hubble rate and p 7 and p v are the photon and neutrino densities, respectively. The photon 
energy density is very well determined from the CMB temperature, so constraints on the expansion in the early 
Universe, H{z), can be translated directly to p v and hence A^ff. 
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Figure 1: A selection of cosmological JV e ff measurements and 68% confidence intervals from the liter- 
ature for various combinations of models and data sets. W denotes WMAP followed by data release. 
The m odels are all ACDM plus th e exten sions given o n the plot. 
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Photon re-heating (electron recombination) 
(T rec « 0.2 MeV, z ICC w 8.5 x 10 8 ) 

When the electrons and positrons decouple from the thermal plasma, most of them annihilate. The released energy 
re-heats the photon population, but not the bulk of the neutrinos since they are no longer in thermal equilibrium 
with the photons. After annihilation 

T v = (4/ll) 1/3 T 7 . (2) 

Because the decoupling of the neutrinos is not instantaneous, some of the high energy neutrinos are still in thermal 
equilibrium with the photons and will feel the re-heating leading to a slight excess of high energy neutrinos relative 
to a blackbody spectrum. The extra energy per particle is well approximated by increasing the number of particles 
slightly, and is commonly absorbed into the definition of N e g via 

Pv=Negz^T?, (3) 

where p„ refers to the neutrino energy density when they are still relativistic, and the factor (7ir 2 /120)T^ is the 
energy density in one species of relativistic neutrinos. The additional re-heating leads to N c s = 3.046 for N v = 3. 



Matter-radiation equality 

(T cq « 0.8 eV, z cq « 3200) 

Relativistic neutrinos behave like radiation, and consequently change the time of radiation-matter equality 
(jLesgourgues fc Pastoll2006l ') 

ao 1 — f v fim 

where £7 m is the total matter density, Q T is the radiation density, and f v = f2„/r2 m is the neutrino fraction. 
The matter-radiation equality determines when the density perturbations can begin to grow and form structures 
through gravitational collapse. If this is pushed to later times, the peak of the matter power spectrum moves to 
larger scales (because the particle horizon is larger when the structures begin to collapse). 



Recombination makes the Universe neutral 

(T rcc » 0.3 eV, z cq » 1370) 

Once the temperature has dropped sufficiently below the ionisation energy of hydrogen, the electrons will couple 
to the nuclei and the Universe becomes neutral (z roc is here defined as when half the electrons are bound). Since 
the neutrinos have already decoupled a long time ago, they are unaffected by this, but the recombination is the 
beginning of the CMB era. 



Photon decoupling and last scattering 

(T* » T ls » 0.25 eV, z* « z ls w 1100) 

When the temperature of the Universe has dropped to m 0.25 eV, the collision rate between photons and free 
electrons drops below the expansion rate and photons propagate freely. The last scattering between photons and 
electrons is observed today as the CMB. The neutrinos do not care about this, but the CMB is one of the most 
powerful probes available to study the early Universe (see Sec. 13. 2[) . 



Baryon decoupling and the drag epoch 

(Tdrag ~ 0.24 - 0.25 eV, z drag w 1020 - 1060) 

Even after the photon-electron scattering ceases and the mean free path of photons becomes effectively infinite, a 
small number of the photons do still interact with the far less numerous baryons (nt,/n 7 ~ 10 -10 ), prolonging the 
pressure wave in the baryons. This is essentially a difference between when the photons stop noticing the baryons 
(decoupling), and when the baryons stop noticing the photons. This delay is known as the drag epoch, at the end 
of which the density fluctuations that form the seeds of structure are frozen in. This event does not affect the 
neutrinos, but the decoupling allows the baryons to cluster in the dark matter potentials initiating the formation 
of observable structures. 



Becoming non-relativistic 

(T nr w 0.05 eV, z m ~ 200 depending on m„) 

At least two of the neutrino mass eigenstates are massive so their velocities will decay adiabatically with the 
expansion of the Univers e, and eventually they become non-relativistic (m„c 2 > p v c) (Lesgo urgues fc Pastorll2006l : 
Ide Bernardis et alJl2009h 

Znr ^2xl0 3 ^. (5) 
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The last eigenstate is allowed to be massless and if this is the case, it will remain relativistic forever. When 
the neutrinos become non-relativistic they behave like a species of warm/hot dark matter, slowing the growth of 
density fluctuations on scales smaller than their free-streaming lengtlQ (see Sec. 13. 3[) . 



3 Observable effects 

The main effect of dark radiation is to alter the expansion rate and the time of matter-radiation equality, which 
affects nucleosynthesis, the CMB, and structure formation. The relevant observables for nucleosynthesis are 
the deuterium and helium abundances, for the CMB it is the power spectrum of temperature and polarisation 
fluctuations expressed in terms of spherical harmonics, and for structure formation it is the (matter) density power 
spectrum inferred from e.g. galaxy surveys. Fig.[3]illustrates the effects of increasing the neutrino mass (solid lines) 
and number of species (dashed lines) on the CMB and density power spectra. In the following sections we will 
discuss these observable effects in detail. 

What is probed directly by cosmological analyses is no t the neutrino mass, but the neutrino density, p v , which 
can be expressed in terms of the sum of neutrino masses (Lcsgourg ues fc Past or 2006) 

" ~ Pc ~ 93.14 eVft 2 ' W 

where p c is the critical energy density for a flat Universe. Because of the small-ness of the measured mass differences 
relative to the upper limits, it is reasonable to assume that if the individual neutrino masses are near the upper 
limit, they are effectively equal. If the individual masses are close to the lower limit, the hierarchy will play an 
important role (jjimenez et alJl201Ch . 



3.1 Big Bang Nucleosynthesis 

The presence of dark radiation mainly alters the nucleosynth esis by increasing the expans i on rate leading t o higher 
deuterium and helium abundances than are observed today (Lcsgourgues & Pastor 2006; Steigman 2012). 



3.1.1 Helium 

The freeze-out of weak interactions determines the neutron/proton ratio, that will eventually determine the helium- 
fraction (Yp = pHe/pbaryon), since all available neutrons either become bound in helium or decay. Increasing the 
expansion rate will induce earlier freeze-out of the weak interaction (see Fig. [3]), and more neutrons will be available 
leading to a higher helium fraction. 

Helium-4 ( 4 He) is the nucleosynthesis relic most sensitive to N e ff but because it is produced in stars it is not 
trivial to obtain a robust value of the primordial abundance. High quoted values of N a g from nucleosynthesis 
are usually driven by the adopted value of Y p an d its uncertainty. Since the beginning of the 1990s the observed 
value has increased fro m Y p = 0.2462 ± 0.0015 l|lzotov et al.lll999h to Y p = 0.2565 ± 0.001(stat.) ± 0.005(syst.) 
(jlzotov fc T huan 2010) which drives the value of iV e ff up. 



3.1.2 Deuterium 

The deuterium abundance is determined by how much deuterium is fused into 4 He after deuterium formation has 
ceased because of too low neutron and proton densities towards the end of the nucleosynthesis. Increasing iVeff 
willl speed up the expansion (but not as much as when helium formed ) allowing for less deut erium to be fused 
into helium and consequently leading to a higher deuterium abundance (Nollctt & Holder 2 0111 '). 



3.1.3 Lithium 

The lithium abundance predicted from standard nucleosynthesis is a factor of three higher than that observed in 
metal poor halo stars. Increasing N e g does not change this problem (Jsteigman 2oT2). 



3.1.4 Results from Nucleosynthesis 

There have been a number of recent attempts to derive N e g from nucleosynthesis alone e.g. Mangan o fc Serpicd 
(|201ll ): lPettini fc Cookd (l2012h . which all seem to be consistent with N e e = 3. Analy ses combining nucleosynthesis 
and CMB appear to prefer N > 3 (e.g. INollett fc Holderffeoill ; iHamann et al.ll201lf h 



^The free-streaming length is the average distance the neutrinos propagate freely, and for kinematic reasons neutrinos 
cannot be confined (or kept out of) regions smaller than the free-streaming length. 
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Figure 2: The redshifts of important events in the early Universe (described in the main text) and their 
dependence on N e g (the first two columns) and on cosmological parameters for fixed N g=A (the second 
and third columns). Changing N c g mainly changes the redshift of matter-radiation equality (z eq ) for fixed 
cosmological parameters. However, this change can be compensated by increasing the values of Hq and 
^m,o- Notice that the neutrino decoupling remains early allowing one to distinguish between the first and 
third scenarios through deuterium and helium abundances. 
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Figure 3: Illustration of the CMB (left) and matter (right) power spectra for a fiducial cosmology and 
how they change for varying neutrino mass (solid lines, we vary /„ = f2„/f2 m = Y, m v /(93. 14 eV£l m h 2 )), 
and effective number of neutrinos (dashed lines) fixing all other parameters. The bottom two plots show 
the ratio of the power spectra when varying Y m v and N e g relative to ^ m v =0 eV, N e g=0. Y m v does 
not affect the CMB power spectrum much, but changes the matter power spectrum. The effect of N e g 
is clearly visible for small scales (high values of I) in the CMB power spectrum, and the two parameters 
are clearly degenerate for the matter power spectrum unless the peak position is very precisely measured. 
The horizontal lines indicate the coverage by current and near future experiments, and the shaded areas 
indicates the relative magnitude of the uncertainties for WMAP (light grey) and SPT (dark grey) in the 
CMB plot, and WiggleZ (dark grey) for galaxy surveys. 
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3.2 Cosmic Microwave Background 

The distribution of hot and cold spots in the CMB is given by the fluctuations in the photon-baryon plasma at 
the surface of last-scattering. These fluctuations are generated by sound waves in the plasma, where the repulsive 
electrostatic force is balanced by the gravitational attraction of dark matter potentials that have been growing 
since z eq (see Fig. [2|. They manifest in the observed CMB power spectrum as actual oscillations, as shown in 
Fig. [3] Accurate measurement of the positions and heights of the peaks and troughs can give us two important 
pieces of information about the relativistic energy density. 

Firstly, the amplitude of the oscillations relates to the matter-to-radiation density ratio. Pressure support 
of the compression wave in the plasma causes the gravitational potential to decay, enhancing the amplitude of 
these perturbations. This enhancement occurs on ly for perturbations that have already crossed the horizon before 
matter-radiation equality l|Hu fc Sugivamalll99a ). which is true for the third and subsequent peaks. This is why 
the third peak is enhanced with respect to the second, and the amplitude ratio of the first to third peaks gives a 
good measurement of z eq , and so N c g. 

Secondly, the peak positions give the size of the sound horizon at decoupling, namely the maximum distance 
the acoustic oscillations can travel before the sound speed in the fluid drops to zero. This relates directly to the 
Hubble rate, and so is also influenced by N e a . Notice from Fig. [3] that increasing the value of N e B moves the 
acoustic peaks to smaller scales (since the sound horizon, r s , becomes smaller for increasing N e g), whereas in the 
matter power spectrum the turnover moves to larger scales (related to the change of z cq seen in Fig. [2]). This 
anti-correlation between effects in the acoustic peaks in the CMB and the turnover in the matter power spectrum 
provides an important consistency check between the early (pre-CMB) and late time (post-CMB) physics. 

In addition, we must consider the damping of the CMB power spectrum 'tail' on small scales. Photons in 
the photon-baryon fluid will diffuse out of over-dense regions, and generate an anisotropic stress that damps the 
oscillations on small scales (i > 200). This effect is called Silk or diffusion damping. The neutrinos contribute to 
the small-scale damping both indirectly - by changing the expansion and thereby the damping from the photons 
- and directly - through their own free-streaming. 

Increasing the expansion rate by adding extra relativistic material will increase the Silk damping because 
structures have less time to grow and are consequently easier to wipe out through diffusion. Similarly, increasing 
the expansion rate increases the helium abundance, which also increases the damping because the helium abundance 
controls the number of free electrons per unit baryon mass. 

Neutrinos contribute to the damping through anisotropic stress and free-streaming effects (|Hu et alj|l995t ). 
Similar to the photons, the neutrinos have a viscosity that creates an anisotropic stress and consequently generates 
damping on small scales. Fluctuations in the neutrino density lead to a shift in the positions of the acoustic peaks. 
These effects are non-negligible, as at the time of decoupling, the standard neutrinos contribute around 10% of the 
total energy density of the Universe. Because the neutrinos are free-streaming, they suppress structure formation 
on scales smaller than their free-streaming length (see Sec. I3.3f) . The combination of these effects tell us something 
more than simply the density of the relativistic species, and can provide details on both the mass and viscosity of 
the dark radiation. 

3.2.1 Results from CMB 

CMB observations of the first three acoustic peaks in the temperature po wer spectra from WM AP have accurately 
measured the redshift of matter-radiation equality to be z cq — 3145^^ (Komatsu ct al. 2009), for models where 
N e H is allowed to vary. However, WMAP cannot measure the tail of the CMB temperature power spectrum 
[I > 1200), and so is unable to put an upper limit on N e $, giving only a lower limit of jV c ff > 2.7 (95% confidence 
level). An upper limit can be found either by combining WMAP with late-time measurements of the matter 
density (SN, BAO, or galaxy power spectra), or through measurements of the CMB tail. When combined with 
small-scale CMB experiments such as the Atacama Cosmology Telescope (ATC) or So uth Pole Telescope ( SPT), 
the constraints improve (see Fig. [TJ, gi ving N eB = 5.3 ± 1.3 from WMAP + ACT (jDunklev et all 1201 if ) and 
N eS = 3.85 ± 0.62 from WMAP + SPT (jKeisler et al.ll201lD . 

In terms of the neutrino mass, the CMB alone cannot provide very strong constraints. In order to have a 
measurable effect on the CMB, the neutrinos would have to be non-relativistic before photon decoupling, but 
current cosmologic al limits indicate that they were relativistic at this time. The upper limit imposed by WMAP 
is ^2 m„ < 1.3 eV (|Komatsu et alj|2009l ). which is close to (in fact beyond!) the absolute limit of ^ m„ < 1.5 eV 
one can theoretically expect to derive from CMB alone (Ichikawa ct al. 2005). However, the CMB does play 
an important role of providing good 'base- level' constraints on the matter density (f2 m ) and the amplitude of 
primordial fluctuations (A^j), which need to be jointly measured when combining the CMB with other cosmological 
datasets to measure the neutrino mass. 

3.3 Large Scale Structure 

Like the CMB, the large scale structure is affected by neutrinos through the expansion and their free-streaming 
effects. The expansion determines the time available for sound waves to propagate in the baryon-photon fluid 
before recombination, which will imprint on the BAO scale measurable from both the CMB and in the galaxy 
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power spectrum (see Fig. [2]). Increasing N e s pushes the BAO peak towards smaller scales, indicated by a slight 
shift of the positions of the 'wiggles' in the galaxy power spectrum in Fig. [3] The position of the matter power 
spectrum turnover is determined by the time available for structures to form, and is consequently pushed to smaller 
scales as matter-radiation equality moves to later times with increasing N e g. However, this effect is degenerate 
with effects from the ordering of the neutrino masses (the hierarchy) that dictates whether one third, or two thirds 
of the neutrinos will b e 'heavy' or 'light'. The light ones will b e radiation-like for longer and enhance structure 
growth on large scales (jjimenez et al-lboiOl ; IWagner et al]|2012fl . 

Massive neutrinos (or other particles) will also slow down structure formation due to their free- streaming prop- 
erties . After decoupling, their masses will lead to an adiabatic decay of their thermal velocities as (jKomatsu et al.l 

noH) 

"thermal = 151(1 + z)(l eV/ra„) kms" 1 . (7) 

Eventually they become non-relativistic and behave as a species of warm/hot dark matter, suppressing growth 
of density fl uctuations on scales charact erised by their free-streaming length at the time when they become non- 
relativistic (Les gourgues fc Pas tor 2006) 



h " = \/i« thcrm f(S)(i + .) (8) 



(1 + z) 2 



This happens because the neutrinos cannot cluster on scales smaller than their free-streaming length and conse- 
quently the density fluctuations are suppressed by a constant factor proportional to the neutrino density, and will 
grow more slowly than without the presence of massive neutrinos. The effect is observable in the density power 
spectrum as a suppression of structure on small scales shown in the right part of Fig. [3] However, the presence 
of neutrinos also enhances the density perturbation s on large scales (by adding to the radiation density for fixed 
f2m,o), which boosts structure formation in general l|Hou et al. 201l]). 

The evolution of the density perturbations can be predicted by solving the coupled set of Einstein and Boltzman 
equations. For small density perturbations the linearised equations are sufficient (calculated numerically e.g. using 
CAME0 or CLASS0), but when the density perturbations grow sufficiently large, the higher-order, non-linear 
terms become important and have to be included, which is computationall y very intensive (and not yet solved 
analytically). Inste ad one has to rely on second order pertu rbation theory (e.g. lSaito et al. 2009; T aruva et al.ll2012h 
or simulations (e.g. lSmith et al . 2003; Je nnings et al.|[2010f ). Num erical solutions show that for /„ = Q, v /Q, m < 0.07 
the suppression is 8P/P = —8f v fo r linear structure formation |H u et aLl ll998T) and the effect increases for non- 
linear structure formation (8 -> 9.6 iBrandbyge et al.ll2008l ; iBrandbyge fe Hannestadll2009l , 120101 ; IViel et al.ll201(il ; 
lAgarwal fc Feldman|[20ll l. 

The neutrinos affect the total matter distribution which is dominated by the dark matter, and consequently 
not directly observable. Instead we use tracers such as galaxies or the absorption lines from neutral hydrogen 
clouds along the line of sight to distant quasars (Lyman alpha forest), for which the power spectrum may differ 
slightly from the total matter power spectrum if baryons do not perfectly trace the dark matte r. Depending on 
the t racer, the bias, b defined as P ODS = b 2 P m , might be scale dependent or scale-independent ([Schulz fc White] 
2006). A scale dependent bias can introduce a degeneracy with the slope change induced by the neutrinos. 

Fo r a mo re in-depth, technical account of massive neutrinos in cosmology, see iLesgourgues fc Pastor! (|2006T ) or 
IWond (|201ll ,. 



3.3.1 Results from large scale structure 

Large scale structure data alone cannot constrain all the cosmological parameters, so they are usually analysed 
j ointly with CMB data. Currently there is no precise measurement of the matter power spectrum turnover 
(jPoole et al.l I2012T) . so large scale structure data mainly improves on the ^ m„ constraints, and are conse- 
quently most important when fitting for ^ m v and N e g simultaneously rather than for N e g alone. The cur- 
rent limits of Y,m v < 0.3 e V for ACDM + Y^ra v (jMoresco et all 120121 : IXia et all [20121: Ide Putter et al] 120121: 
Riemer-Sorensen et al.ll2012r) ^ and Y m v < 0.65 eV fo r ACDM + J2m v + N B « fe.g. IHamann et al.ll2010l : Ijoudakil 
20121 ; iRiemer-Sorensen et al]|2013l ; IWang et alj|2012T ) are stronger than what is provided by particle physics ex- 
periments. 



s http : //camb. inf o| 

•- http : //lesgourg . web . cern . ch/lesgourg/class . 
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Figure 4: The 68% and 95% confidence limit contours for fitting CMB (WMAP7+SPT) + large scale 
structure (WiggleZ) + H[z) + BAO (SDSS, 6dFGS, BOSS) + SN (SNLS) data with ACDM (dashed red) 
and extensions to t he model with m v (dashed blue), N e g (dot-dashed green), and both m„ and N e g 
(solid black) from (jRiemer-Sorensen et al.ll2013lh The resulting contours are consistent with each other 
for all parameters clearly showing that the preference for high values of N c g is not related to an off-set in 
the measurement of a single standard ACDM parameter, but can be added within the standard parameter 
uncertainties. 



4 Discussion 



4.1 Do we need N eS ? 

Perfect measurements of the standard ACDM parameters do not exist. Adding N c g to the cosmological model 
does not shift the best fit values of the standard ACDM parameters in any systematic way, indicating that the 
preference for high N e s is not due to an off-set measurement of a single parameter, but can be accommodated 
within the uncertainties of the standard parameter values. 

Fig. 2] compares the contours from fitting a range of cosmological data (see caption) with a standard ACDM 
model alone, and with yj m„ and N e B extensions individually and combined. Naturally the constraints are tighter 
for pure ACDM, since there are fewer parameters to constrain, but overall the constraints are consistent for all 
parameters: Adding yj m v to the model (dotted blue) does not change the ACDM parameters or uncertainties 
significantly whereas adding N e s (dot-dashed green) increases the values and uncertainties, particularly for Ho- 
Adding both parameters simultaneously (solid black) increases the preferred values of Ho and n s but also the 
uncertainties, so the values remain consistent with the pure ACDM case. The shift in n s happens because the 
increased expansion (Ho) changes the height of the peaks in the CMB, as it corresponds to increasing the physical 
matter density. Increasing N e « recovers the details of the CMB peaks (because the original ratio of matter 
to radiation is recovered, restoring the details of the Silk damping), but with too much power on large scales, 
which can be eliminated shifting the primordial power spectrum from very red (n s = 0.96) to slightly less red 
(n s = 0.98 - 1.0). 
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4.2 Priors 

Since most constraints on N e g are obtained through Bayesian analyses, one co uld be concerned that the prefe rence 
for N e ff> 3 orig inates in the choice of priors on the model parameters. Both [Gonzalez- Morales et al.l l|201ll ) and 
iHamannl l|2012f ) investigated this and found no significant shifts in the N c g constraints for various selections of 
priors. 

4.3 Beyond ACDM+E^+iV cfr 

Other physical effects such as curvature, varying equation of state, runnin g of the spectral index etc. could mimic 
N c g> 3 if not properly accounted for in the modelling. Ijoudakil (|2012h demonstrated that the deviation from 
N e ff — 3.046 is diminished if allowing for curvature, varying equation of state, running of the spectral index, 
and/or the helium fraction by fitting to a combination of CMB + BAO + HST + SN data. However, for all 
parameter combinations the preferred value of N e g was still above three (see Fig. [T]), and only when more than 
one extra parameter extension of the ACDM cosmology (e.g. curvature and varying equation of state in addition 
to Ngfs and ^2 m v ) were considered, the preferred value of N B s became consistent with three within one standard 
deviation. Consequently the preference for N e s =/ 3 does not appear to be strongly d ependent on the complexity 
of the cosmological model. A similar conclusion was reached by iMenegoni et all (|2012h which concluded that even 
when varying both the helium fraction and the fine structure constant on top of ACDM + N c s the preferred value 
of i Veff remained high. 

lAudren et al.l (|2012h took an alternative approach and separated cosmological effects originating directly from 
early Universe conditions from those related to later time effects, when analysing CMB data. Marginalising over 
the late time effects, they found the preferred value of N e s to be consistent with three, and consequently suggested 
that the preference for high values of N c g does not relate to changed expansion in the early Universe, but rather 
to our understanding and modelling of late time effects. 

4.4 Thermalisation 

All current analyses relating the dark radiation to neutrinos rely on an assumption of the sterile neutrinos being in 
thermal equilibrium with the photons before the neutrinos decouple. However, this is an ad-hoc assumption and 
a full thermalisation calculation inv olving momentum-dep endent density matrices for three active and AN sterile 
flavours has never been performed (|Abazaiian et al.ll2012ft . 

The thermalisation can b e suppressed in the presence of large lepton asymmetry (more leptons than anti- 
leptons, iFoot fe Volkasl fl995l ). In that case the mass constraints of y^.m v < 0.3 eV from l arge scale structure 
do not apply because they rely on the assumption of thermalisation ( Hannestad et al"1l2012f ). This allows for a 
joint interpretation of the cosmology and particle physics results under very specific early Universe conditions. 
The necessary lepton asymmetry (L = (n/ — nj)/n 1 « 10~ 2 ) is below any observational constraints, but well 
above the usually assumed value which is of the same order of magnitude as the baryon asymmetry [L ~ 10 -10 ). 
Large n eutrino asymmetr ies can be generated by the active-sterile neutrino mixing without the need to invoke new 
physics (|Foot et al.lll996h . 

4.5 Non-standard decoupling 

The number of relativistic degrees of freedom depends on the decoupling processes involved. In the three-neutrino 
framework, the non-instantaneous decoupling of the standard model neutrinos leads to an expected value of 
^V e ff=3.046 for standard decoupling. However, non-standa rd decoupling can maximally increase this value to 
iV e ff=3.12 for three neutrino species (Mangano ct al. 2006), and consequently non-standard decoupling cannot 
alone explain the preference for the currently observed high values of N e g. 

4.6 Particle properties 

The value of N e a> 3 is often interpreted as an extra specie of non-interacting relativistic neutrino-like particle. 
If this is the case, it should behave as a neutrino during the evolution of density perturbations which is can 
be characterised b y a relation between its effective sound speed and viscosity parameter of c\^ = c^ is = 1/3 
(jSmith et al.ll2012l ). c v i s describes the ability of the particles to free-stream out of gravitational potentials, and 
Ccff 7^ 1/3 allows for the existence non-adiabatic pressure perturbations. If the dark radiati on is composed of a 
particle with significant interaction, these parameters can take other values l|Basbe>ll et alj2009fl. thus allowing us to 
distinguish between hypothesised species (|Trotta fc Melchiorrill2005h . This was tested bv lArchidiacono et alj (|201lf ) 
using a combination of CMB + large scale structure + HST data. They found N c s > 3 at 2o (N e g = 4.I5I023) and 
c cff — 0.241q'q2 consistent with 1/3 while c^ is remains unconstrained. This result is consistent with a relativistic 
non-interacting particle but is not a unique signature of such a particle. 

For the specific case of hypothetical sterile neutrinos there is another effect to consider, which could change 
the nucleosynthesis, namely the distortion of the electron-neutrino phase space distribution through neutrino 
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oscillations. For the parameters favoured by short baseline experiments this effect is negligible because the sterile 
and active neutrinos will be in thermal equilibrium and consequently share the same phase space distribution. 

4.7 Time evolution of N e g 

There have been a nu mber of attempts to derive N e ff from nucleosynthesis alone (e.g. Mangano fc Serpicoll20lil ; 
iPettini fc Cookd 12012? ) that all seem t o be consistent with N e g = 3, while analyse s combining nucleosynthesis 
and CMB seem to prefer N > 3 (e.g. iNollett fc Holderl I20TH ; lHamann et all 120111 ). It is concerning that the 
preference for N e « > 3 is present in all analyses including CMB, but not preferred by nucleosynthesis alone. This 
could indicate a systematic error in one of the data sets. However, nucleosynthesis alone relies on very few and 
notoriously difficult measurements of the deuterium and helium fractions. 

The difference between the results could also be interpreted as potential tension between N e g measured at 
two different epochs (nucleosynthesis and recombination), where we naively would expect N e a to be the same. 
A temporal variation could be explaine d theoretically by particles decaying to photons or neutrinos and thereby 
altering their temperature distributions (Ichikawa ct al. 2007; Fischlcr & Meyers 2011; Eggcrs Biacldc 'et al.ll2012l ; 
iHooper et alll2012l : iMenestrina fc Scherrerll2012f l. or heavier particles in thermal equilibrium with the neutrinos, 
but not the photons l|Boehm et al,ll2012h . Currently there is no experimental evidence for the existence of such 
particles. 

4.8 Link between particle physics and cosmology 

Is the cosmological dark radiation related to the sterile neutrinos from the oscillation experiments? 

The neutrino oscillation results favour a mass of order 1 eV, which is incompatible with cosmological mass 
constraints from combinatio ns of CMB and large scale structure measurements (e.g. ^ m„ < 0.3 — 0.6 eV, 
lRiemer-S0rensen et al.l [20121 ). lArc hidiacon o et all (|2012l ) analysed the short baseline oscillation data together 
with cosmological data and found that short baseline data and CMB are compatible with both 3+1 and 3+2 
scenarios, but when adding large scale structure both scenarios are under serious pressure. However, it has been 
shown that factors such as initial lepton asymmetry (Sec. 14. 4p can alleviate these constraints by introducing a 
non-thermal neutrino s pectrum for which the neutrino mass constraints from large scale structure can be avoided 
(jHannestad et al.ll2012T ). 



5 Future measurements of N e g 

Due to its high spatial resolution and several spec tral bands, the Planck microwave observatory, is expected to 
meas ure N a g with a ler uncertainty of 0.20 — 0.25 (B ashinsky & Scliak 2004; Hanncst ad et al.ll2006l ; iPoole et all 
I2012T ). The neutrino anisotropy will also lead to a polarisation of the CMB detectable in future experim ents. Com- 
bining Planck with future small-scale polar isation experiments, a sensitivity of 0.05 may be achievable (|Galli et all 
l2010l ; lThe COrE Collaboration et alj|201 ll ). which would almost allow for a measurement of the thermal distortion 
of the neutrino spectrum (the 0.046). 

Measuring the posit ion of the peak of the matter power spectrum (the turnover) would give another handle on 
N e f[. IPoole et al.l l|2012h predicts that a Euclid- like galaxy survey will be able to constrain N e g with a precision of 
20% in dependently of the C MB. However, the position of the turnover can be degenerate with neutrino hierarchy 
effects l|Waener et al.ll2012h . The degeneracies between ^m„, A'eff, and the hierarchy allows for a measurement 
of either of them if the remaining parameters can be 'fixed' by independent data. 

Any neutrino-like behaving particle (including sterile neutrinos and axions etc., which decouple early as rel- 
ativistic and become non-relativistic) can mimic the effect of the neutrinos. If the value of N e g > 3 is due to a 
particle, the result points towards physics beyond the Standard Model, and its existence it will have to be con- 
firmed by particle physics laboratory experiments. If no compatible particles are found in laboratory experiments, 
the cosmological preference for additional species indicates a lack of understanding and modelling of cosmological 
data, both of which are interesting scenarios. 



6 Conclusions 

Cosmological data seems to favour the presence of dark radiation at the time of photon decoupling at a 2cr level, 
while short baseline particle physics experiments are best explained by the existence of 1-2 sterile neutrinos. 
Unfortunately the details and required particle properties of the two scenarios are hard to reconcile, and appear 
to have separate solutions. No matter whether the cosmological trend is due to new particles or systematics, it is 
a clear indication of something we do not understand. 
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